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Abstract - If supersymmetry (SUSY) has a compressed spectrum then the current mass limits 
from the LHC can be drastically reduced. We consider a possible 'worst case' scenario where the 
gluino and/or squarks are degenerate with the lightest SUSY particle (LSP). The most sensitive 
searches for these compressed spectra are via the final state LSPs recoiling against initial state 
radiation (ISR). Therefore it is vital that the ISR is understood and possible uncertainties in the 
predictions are evaluated. We use both MLM (with Pythia 6) and CKKW- L (with Pythia 8) 
matching and vary matching scales and parton shower properties to accurately determine the 
theoretical uncertainties in the kinematic distributions. All current LHC SUSY and monojet 
analyses are employed and we find the most constraining limits come from the CMS Razor 
and CMS monojet searches. For a scenario of squarks degenerate with the LSP and decoupled 
gluinos we find Mq > 340 GeV. For gluinos degenerate with the LSP and decoupled squarks, 
Mg > 500 GeV. For equal mass squarks and gluinos degenerate with the LSP, Mq i3 > 650 GeV. 
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Introduction. — Supersymmetry (SUSY) 1 has now 
been extensively searched for at the LHC [2148]. In the 
limit that the lightest supersymmetric particle is mass- 
less the bounds on equal mass squarks and gluinos are 
> 1.5 TeV. However, if the LSP is massive and degener- 
ate with the squarks and/or gluinos, these limits can be 
severely weakened [9||l3] . 

In the case of 'extreme' compression, initial state radia- 
tion (ISR) in the form of jets may be used to discover the 
model [14]. This method has now been used to set mass 



limits on squarks and gluinos Illlll2] and stops 15 16 



Other ideas to look for compressed SUSY include using 
monophotons 17 and soft leptons Il8]. Here we look 
at simplified models in which the first two generations of 
squarks and/or the gluinos can be quasi- degenerate with 
the LSP and set limits on the particle masses. These mod- 
els represent a 'worst case' scenario for R-parity conserv- 
ing SUSY at the LHC and thus give a model independent 
bound on these masses. 

The limits completely rely on the accuracy of the ISR 
prediction for radiated jets. Thus, we use two independent 
matching schemes (CKKW-L and (MLM [2l]) to 



carefully understand the accuracy of the prediction. In 
addition, matching scales and parton shower variables are 
changed to ensure that the limits are robust. 

ISR and matching procedure. — To reliably set 
limits in SUSY models with compressed spectra it is vital 
that the ISR is well modeled and any uncertainties are 
understood. To achieve this we must use a matrix element 
prediction for hard and well separated partons whilst using 
a parton shower for soft and/or collinear jets. To bring the 
two approaches together we must match the two methods 
in a consistent algorithm. 

The reason that a matrix element prediction for hard 
radiation is required is that a parton shower approach con- 
tains large uncertainties in this regime. They are mainly 
due to the choice of starting scale for parton shower evolu- 
tion. This has been shown to have a large effect in squark 
and gluino production [22] and we show the variation be- 
tween the 'wimpy' and 'power' shower settings in Pythia 
6 [23] and 8 [24] in Fig[l] 

In order to match the matrix element to the parton 
shower we must be careful to avoid double counting so 
that areas of phase space are only filled by one approach. 
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Fig. 1: A comparison of the uncertainty in the ISR jet pr distribution for the production of squarks (Mq = 500 GeV) between 
the parton shower prediction (green, light), MLM matching (pink, medium) and CKKW matching (blue, dark). The parton 
shower uncertainty is found by varying the Pythia 6 and 8 parton showers between the 'wimpy' and 'power' settings [22] . 
The matching uncertainties are found in both cases by varying the matching scales between 50 and 200 GeV and additionally 
for MLM matching by varying the parton shower between the 'wimpy' and 'power' settings. For the softer jet 3 (the first 
unmatched jet), the relative uncertainty of the parton shower approach is reduced since the phase space for this emission is 
better constrained. 



In addition we would like a smooth transfer between the 
different areas of validity. Finally, the prediction should 
not have a significant dependence on the chosen match- 
ing scale or parton shower. Within SUSY we have the 
additional problem that we can double count events with 



resonant propagators, which must be removed 22,25 



To independently check the predictions of the match- 
ing algorithm, two approaches were used. First is the 
integrated MLM [2l] matching available in Madgraph 5 
[25][26] which is interfaced with the Pythia 6 shower [23] . 
To estimate the uncertainty, we varied the matching scale 
between 50 and 200 GeV and independently the parton 
shower between the 'wimpy' and 'power' settings. As seen 
in FigJT] this results in a significant reduction in the uncer- 
tainty compared to the parton shower alone. The second 
approach was the CKKW-L 19 , 20] matching algorithm, 



developed for Pythia 8 24,27 . It was adapted to work 
with SUS^Qand gives consistent results to those obtained 
with MLM matching, Fig[l] 

Simplified Models. — In order to reduce the SUSY 
parameter space and place model independent limits, we 
use three simplified models. The idea is to investigate 
the 'worst case' scenario for R-parity conserving SUSY. 
We thus assume that either the first two generations of 
squarks or the gluinos or both are quasi- degenerate with 
the LSP. The degeneracy has the effect of making all of 
the SUSY decays invisible to the detector as the produced 
charged particles are too soft to be reconstructed. There- 
fore, events with ISR are solely relied upon to set any 
limits on the model. 



Our first scenario is labeled the 'Decoupled Gluino' 
model, Fig[2|a). Here the first two generations of squarks 
are quasi-degenerate with the LSP (1-100 GeV mass split- 
ting) while the gluino is completely removed from the sce- 
nario. The idea is to set a lower limit on the first two 
generation squarks masses that is completely independent 
of the gluino mass. The third generation of squarks are left 
free (obviously heavier than the LSP) because in general 
the Yukawa contribution to the running of the mass leads 
to a splitting between these and the first two generations 
of squarks. However, a degenerate contribution can easily 
be added by simply rescaling the cross-section by 5/4 for 
only sbottoms or 6/4 for stops as well. 

The second scenario we name the 'Decoupled Squark' 
model, Fig[2|b). Here the gluino is quasi-degenerate with 
the LSP (1-100 GeV mass splitting) and the first two gen- 
erations of squarks are removed from the model. In the 
limit that all squarks are removed from the scenario it 
must be stated that the gluino becomes stable and a dis- 
tinctive signal would therefore be seen as so called 'R- 
hadrons'. In fact, even for moderate squark masses it is 
easy to make a gluino in a compressed spectra long-lived. 
However, it is possible that the third generation squarks 
could be much lighter than the other squarks. These could 
mediate prompt gluino decay whilst having a negligible 
impact on the search. Therefore we assume a prompt de- 
caying gluino in this scenario as an interesting limiting 
case0 

As a third scenario we consider the 'Equal Mass' model, 
Fig. |2|c). Here, the gluino mass is set quasi-degenerate 
with the LSP (1-100 GeV mass splitting) and the first two 



! We would especially like to thank Stefan Prestel for his invalu- 
able help in adapting the algorithm. 



2 Such mo dels are already being investigated by the LHC collab- 
orations. |28l 
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Fig. 2: The spectra for the simplified models studied in this paper. For the 'Decoupled Gluino' scenario we remove the gluino 
from the spectrum and vary the mass difference between the squarks and the LSP from 1 to 100 GeV. In the 'Decoupled Squark' 
scenario we remove any squarks from the spectrum and vary the mass difference between the gluino and the LSP from 1 to 
100 GeV (see accompanying text for caveats). In the equal mass scenario, the gluino is set as the most massive particle and 
the squark mass is halfway between the gluino mass and LSP mass. In all scenarios the squarks are a summation over the first 
two generations and the third generation are ignored. 



generations of squarks have a mass halfway between the 
gluino and the LSP, 



Mq = ^(M~ g + M LSP ). 



(1) 



The third generation squarks are once again ignored but 
there is no conceptual reason why they could not be added. 
However, they would only give a very small contribution to 
the final cross-section so our limit would remain practically 
unchanged. 

Additional assumptions have to be made on the mod- 
els in order to use the LHC analyses. Firstly, we assume 
that all decays are prompt and we have no new heavy 
states traversing the detector or producing displaced ver- 
tices. In any case, these signals are currently searched for 
and should provide a distinctive signature [29] . Secondly, 
when the mass splittings between the squarks/gluinos and 
the LSP is increased, we assume that the decay occurs 
directly to the next lightest particle in our decay chain. 
Therefore, we assume that no other states exist in be- 
tween. In the limit of degeneracy, even if new states did 
exist, the phenomenology is unchanged because all the 
momentum of the initial parent particle is still transfered 
to the LSP. For increased splitting it is possible that the 
limits may be changed with multiple decays to many soft 
particles. However, we do not consider this possibility. 

Searches. — In order to set the most stringent bound 
on the model parameter space we use all current SUSY 
hadronic searches from both ATLAS and CMS. In addition 
we apply the mono jet searches from both experiments. 

The motivation for using all searches is that there exists 
a range of strategies to look for SUSY at the LHC and it 
was not obvious which of these would be most productive 
for compressed spectra. The 'vanilla' LHC searches are 
based around an effective mass and missing energy cut. 
Various signal regions are defined with different amounts 
and proportions of effective mass (or Ht) and missing en- 
ergy. Both ATLAS [2] and CMS [6] have an analysis of 



this kind and they give similar results for mSugra. In ad- 
dition, CMS has a number of 'shape' based analyses using 
a T [5)[30], M T2 [7)[3l] and Razor §[32]. 

As ISR is vital for discovering SUSY in compressed spec- 
tra we would expect that the signal is dominated by the 
hard emission of a single parton from the initial state. 
Therefore, we also use the monojet analyses (ATLAS 2 
and CMS 3 ) that are optimised for precisely this kind of 
signal. 

All the searches were implemented within the RIVET 
[33] analysis package. As we are only interested in jets, no 
experimental efficiencies apart from quality cuts were re- 
quired. Momentum smearing of the jets, including adding 
mis- measured tails were included [34^] However, these ef- 
fects were found to have a negligible impact on all of the 
search efficiencies. All analyses were tested against any 
mSugra or simplified model bounds found in the original 
studies. In addition, where distributions and cut flows are 
available, these were also tested. Acceptances were re- 
produced to within 20% for all analyses but were usually 
found to agree much better. 

In order to set limits in a consistent way across all anal- 
yses, we use the 95% confidence limit given by the Rolke 



Test 35 - 37 for the most constraining signal region in each 



analysis. To be able to fairly compare the analyses and 
gain a better understanding of the optimum search re- 
gions, if the limit is better than the limit expected to have 
been found, we take the expected limit. In this way we 
find more conservative limits than the official analyses. 

Both CMS olt and Razor searches use a combined like- 
lihood over all signal regions to set a limit. We instead 
conservatively use only the best signal region as this allows 
a fairer comparison with the other searches. In addition, 
CMS Razor uses an unbinned likelihood that is impossible 
to reproduce as the data is not publicly available, however 
a 60 bin data set is available. We combined adjacent bins 



3 We thank the authors for providing the details of the jet mis- 
measurement constants in a private communication. 
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Fig. 3: Ratios of the signal cross-section (asig) to the cross-section required for exclusion (vexc) at the 95% confidence level for 
squark masses in the decoupled gluino scenario, Fig|2ja). The monojet searches (left) and SUSY searches (right) are shown. 
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Fig. 4: Limits at the 95% confidence level for squark masses in the decoupled gluino scenario, Fig[2ja), for both monojet (left) 
and SUSY searches (right). The mass splitting between the squark and the LSP is varied between 1 and 100 GeV. 



in this set to give a final search with 20 signal regions 
which is conservative compared to the official analysis. 

At the time of release both the ATLAS monojet and 
CMS aT only had analyses with 1 fb -1 and 1.14 fb -1 . In 
order to give a fair comparison with the other analyses 
we extrapolated these searches to 5 fb -1 by reducing sta- 
tistical errors. This is a conservative approach as many 
systematic errors are also statistically limited and likely 
to be improved. 

Limits. — To calculate limits in our simplified mod- 
els, cross sections are calculated to next-to-leading order 
including known next-to-leading-logarithmic corrections 
using NLL-Fast 38] [40 



The theoretical uncertainty is 
calculated including the factorisation and renormalisation 
scale and PDF 41 errors. In addition, we vary the match- 
ing scale and parton showers, Fig{l] and take the result 
with the least constraining bound. 



We find that in the decoupled gluino scenario with a 
quasi-degenerate LSP, Figj2|a), we can set a bound on 
the first two generations of squarks of Mq > 340 GeV. In 
this scenario we find that the best limit is set by both the 
CMS monojet and CMS Razor searches, Fig(3j whilst the 
other SUSY searches give slightly weaker limits of Mq > 
300 GeV. All the limits along with the relevant search 
regions are given in Tab [I] 

Interestingly, we also see that the extrapolated ATLAS 
monojet search provides a noticeably worse limit of Mq > 
270 GeV. Our analysis showed that the difference between 
the two monojet searches is primarily down to the second 
jet veto used in the ATLAS search which reduces the signal 
acceptance. 

For all the SUSY analyses, we find the best search re- 
gions are where the proportion of missing energy compared 
to the total energy in the event is the highest. This is to 
be expected as the relevant events in our models that are 
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Fig. 5: Limits at the 95% confidence level for gluino mass in the decoupled squark scenario, Fig|5Jb), for both monojet (left) 
and SUSY searches (right). The mass splitting between the gluino and the LSP is varied between 1 and 100 GeV. 
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Fig. 6: Limits at the 95% confidence level for gluino masses in the equal mass squark-gluino scenario, Fig[2|c), for both monojet 
(left) and SUSY searches (right). The mass splitting between the gluino and the LSP is varied between 1 and 100 GeV whilst 
the squark mass is placed between the two. 



expected to have a topology with a dominant jet that will 
recoil against the LSPs. 

As we increase the mass splitting, Mq — Mlsp, we see 
that the monojet searches rapidly lose their effectiveness, 
Fig. [4j This is due to the fact that both the ATLAS and 
CMS searches have a third jet veto. As the mass splitting 
is increased, the SUSY decays produce additional jets and 
these events are vetoed. 

In contrast, the SUSY searches become more effective 
as the mass splitting is increased. For example, the CMS 
Razor limit improves from Mq > 340 GeV when the mass 
splitting is 1 GeV to Mq > 400 GeV when the mass split- 
ting is 100 GeV. This is because the SUSY searches often 
allow any number of additional jets and include the extra 
radiation in the search variable. 

In the decoupled squark scenario, Fig{2|b), we find that 
we can set a bound on the gluino mass, M~ g > 500 GeV, 



Fig|5j and CMS Razor sets the best limit. Finally, in 
the equal mass squark gluino model, FigJ5|b), we find a 
limit of Mq ~ Mg > 650 GeV which is set both by CMS 
Razor and the CMS monojet search, Figj6] As before, 
the limit in the monojet search is reduced as we increase 
the mass splitting but the majority of the SUSY limits 
actually improve and for 100 GeV mass splitting, CMS 
Razor gives a limit, Mg > 700 GeV. 

Conclusions. — We have set limits in simplified SUSY 
models with compressed spectra. We find a limit on the 
first two generations of squarks, Mq > 340 GeV, in a 
model with a decoupled gluino. With decoupled squarks 
we find a limit on the gluino mass of, Mg > 500 GeV. 
With equal mass gluinos and squarks we find a limit of 
Mq ~M~ g > 650 GeV. 

We would like to comment that the results for the de- 
coupled squark and equal mass squark-gluino models are 
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Table 1: Comparison of the bounds on the mass of SUSY particles for the different searches employed at the LHC. The 
luminosity of the searches and the most constraining search region are also given (the search region names refer to those given 
in the original experimental papers). *The ATLAS and CMS monojet searches only give these bounds for mass differences 
< 5 GeV. For larger mass splittings, the bounds become much weaker. + The ATLAS monojet and CMS olt searches are both 
extrapolated from ~1 fb~ x to give a more direct comparison. 



in good agreement for the ATLAS SUSY search with those 
previously presented in 
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